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A CENTRAL LIMIT THEOREM FOR STOCHASTIC RECURSIVE 
SEQUENCES OF TOPICAL OPERATORS 

By Glenn Merlet 

CEREMADE, Universite Paris-Dauphine 

Let (A n )n£N be a stationary sequence of topical (i.e., isotone 
and additively homogeneous) operators. Let x(n,xo) be denned by 
a:(0, xo) = xq and x(n+l,xo) = A n x(n, xo). It can model a wide range 
of systems including train or queuing networks, job-shop, timed dig- 
ital circuits or parallel processing systems. 

When (A n ) n £fi has the memory loss property, (x(n, xo)) n e.N satis- 
fies a strong law of large numbers. We show that it also satisfies the 
CLT if (A n ,) ne N fulfills the same mixing and integrability assump- 
tions that ensure the CLT for a sum of real variables in the results 
by P. Billingsley and I. Ibragimov. 

1. Model. An operator A : M. d — > R rf is called additively homogeneous if it 
satisfies A(x + al) = Ax + al for all x £M. d and a £ R, where 1 is the vector 
(1, . . . , 1)' in R . It is called isotone if x < y implies Ax < Ay, in which the 
order is the product order on R rf . It is called topical if it is isotone and 
homogeneous. The set of topical operators on R rf will be denoted by Top d . 

We recall that the action of matrices with entries in the semiring R max = 
(R U {— oo}, max, +) on R^ax i s defined by (Ax)i = max. j(Aij + xj). When 
matrix A has no — oo-row, this formula defines a topical operator, also de- 
noted by A. Such operators are called max-plus operators and operators 
composition corresponds to the product of matrices in the max-plus semir- 
ing. 

Let (A n ) n& ^ be a sequence of random topical operators on R rf . A stochas- 
tic recursive sequence (SRS) driven by stochastic recursive sequence is a se- 
quence (X n ) n£ j$ satisfying equation X n+ \ = A n X n . To study such sequences, 
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we define (x(n, xo))neN by 

x(0,x ) = Xq, 

(1) 

x(n + 1, xq) = A n x(n, x ). 

This class of system can model a wide range of situations. A review of 
applications can be found in the last section of [5]. When the x(n,-)'s are 
daters, the isotonicity assumption expresses the causality principle, whereas 
the additive homogeneity expresses the possibility to change the origin of 
time. (See Gunawardena and Keane [18], where topical functions were intro- 
duced.) The max-plus case has, for instance, been applied to model queu- 
ing networks (Mairesse [26], Heidergott [20]), train networks (Heidergott 
and De Vries [21] and Braker [9]) or job-shop (Cohen, Dubois, Quadrat 
and Viot [11]). It also computes the daters of some task resources models 
(Gaubert and Mairesse [16]) and timed Petri nets including events graphs 
(Baccelli [1]) and 1-bounded Petri nets (Gaubert and Mairesse [15]). The 
role of the max operation is to synchronize different events. For develop- 
ments on the max-plus modeling power, see Baccelli, Cohen, Olsder and 
Quadrat [2] or Heidergott, Olsder and van der Woude [22]. 

To clarify things, let us introduce a simple example. 

Example 1.1. Our process assembles two parts. The nth time it is done, 
it takes time 03(71). The parts are prepared separately, which respectively 
takes times a\(n) and a 2 (n). Then, they are sent from the preparation places 
to the assemblage place, which takes times t\(n) and t 2 (n) respectively. Once 
the assembly place has finished an operation, it asks for new parts. At that 
time, if a preparation place has a ready part, it sends it and starts preparing 
another one. Otherwise, it finishes the one it is processing, sends it, and 
immediately starts preparing another one. This is summed up in Figure 1. 
We denote by (X n )i and (X n )2 the starting date of preparation of the nth 
part of each type and by (X n ) 3 the starting date of the (n — l)th assembly. 

Sequence (X n ) n ^ is ruled by equations 

(X n+ i)i = max((X n )i + ai(n), (X n ) 3 ), 
(X n ) 2 = ((X n ) 2 + ai (n),(X n ) 3 ) 

and 

(X n+1 ) 3 = a 3 (n) + max((X n )i + h(n), (X n ) 2 + h{n)) 
= max((X„)i +ii(n) + o 3 (n- 1) +ai(n), 
{X n ) 2 + t 2 (n) + a 3 (n - 1) + a 2 (n), 
(X n ) 3 + ti(n) Vt 2 (n)+a 3 (n- 1)), 
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Fig. 1. A simple production system. 



in which we recognize equation (1), with A n defined by the action in the 
max-plus algebra of 

(a\{n) — oo a-sin) 

— oo Ct2( n ) cis(n) 

ai{n) + t\{n) a 2 (n) + t 2 (n) t\ (n) V ti (n) + a 3 (n — 1) , 

We assume that the sequence (A(n)) n ^ is stationary and ergodic. 



We will focus on the asymptotic behavior of (x(n, -)) n eN- It follows from 
Theorem 2.1, due to Vincent, that V« x i{ n ^°))n£n converges to a limit 7. 

In many cases, if the modeled system is closed, then every sequence of 
coordinate (xi(n, X°)) n ^ also tends to 7, by Theorem 2.2. The so-called 
cycle time 7 is the inverse of the network's throughput or the inverse of the 
production system's output, as in Example 1.1. Therefore, there have been 
many attempts to estimate it (Cohen [12], Gaujal and Jean-Marie [17], Resig 
et al. [30]). Even when the A^s are i.i.d. and take only finitely many values, 
approximating 7 is NP-hard (Blondel, Gaubert and Tsitsiklis [7]). Hong and 
his coauthors have obtained [3, 4, 14] analyticity of 7 as a function of the 
law of Ai. They did so under the so-called memory loss property (MLP) 
introduced by Mairesse to ensure some stability of (x(n, •))neN (see [26]). 

We prove another type of stability under the same assumptions. If (A n ) n£ fq 
has the MLP, then (x(n, -)) n ^ actually satisfies a central limit theorem 
(CLT) under the same mixing and integrability hypotheses as the real vari- 
ables in the CLT for sum of stationary variables by Billingsley [6] and Ibrag- 
imov [24]. 

As far as we know, two CLTs have already been proved for this type 
of sequences: one in [30] and the other in [28]. The most obvious improve- 
ment is that both assumed that the A n 's are i.i.d. Moreover, in [30], the 
hypothesis is difficult to check and the A^s are max-plus operators defined 
by almost surely bounded matrices. In [28] the main hypothesis is also the 
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MLP, but integrability hypotheses were stronger, except for a subclass of 
topical operators. 

The remainder of this article is divided into two sections. In Section 2 
we define the memory loss property, present some law of large number type 
results and state our central limit theorems. In Section 3 we prove the theo- 
rems. First, we state the CLT for subadditive processes by Ishitani [25], then 
we check that (\J i Xi(n, 0)) n eN satisfies each of its hypotheses. To this aim, 
we use Mairesse's construction of the stationary version of the SRS, as well 
as different results from ergodic theory, depending on the hypothesis. We 
eventually deduce the results on (x(n, -))neN from those on Xi(n, 0)) ne jq. 

2. Presentation. 

2.1. Memory loss property. Dealing with homogeneous operators, it is 
natural to introduce the quotient space of M. d by the equivalence relation ~ 
defined by x ~ y if x — y is proportional to 1 = (1, . . . , 1)'. This space will 
be called projective space and denoted by P]R^ ax . Moreover, x will be the 
equivalence class of x. 

The function x h-> (xi — Xj)^ embeds PR^ax onto a subspace of rW** -1 ))/ 2 
with dimension d— 1. The infinity norm of R( d ( rf_1 ))/ 2 therefore induces a dis- 
tance on PR^ ax which will be denoted by 5. A direct computation shows that 
S(x,y) = V \{x,i — yi) + \fi(yi — Xi). By a slight abuse, we will also write S(x,y) 
for 8(x,y). The projective norm of x will be \x\-p = S(x,0) = \J i Xi — Ai x i- 

Let us recall two well-known facts about topical operators. First, a top- 
ical operator is nonexpanding with respect to the infinity norm (Crandall 
and Tartar [13]). Second, the operator it defines from Pl^ to itself is 
nonexpanding for 5 (Mairesse [26]). 

The key property for our proofs is below: 

Definition 2.1 (MLP). 

1. A topical operator A is said to have rank 1 if it defines a constant operator 
on PMmax : Ax does not depend on x € R d . 

2. A sequence (A n ) n ^ of Top d -vedued random variables is said to have the 
memory loss property (MLP) if there exists an N such that An ■ ■ ■ A\ 
has rank 1 with positive probability. 

This notion has been introduced by Mairesse [26], with the ^4 n 's as max- 
plus operators. In this case, the denomination rank 1 is natural. 

We have proved in [27] that this property is generic for i.i.d. max-plus 
operators: it is fulfilled when the support of the law of A\ is not included in 
the union of finitely many affine hyperplanes, and in the discrete case the 
atoms of the probability measure are linearly related. 
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This result applied to Example 1.1 states that the sequence (A n ) n ^ has 
the MLP provided that the support of (ai (2), 02(2), ti (2), £2(2), 03(1)) is not 
included in a union of finitely many affine hyperplanes of IR 5 . This is not 
completely straightforward because the matrix A(l) is defined by only 5 
variables, but the detailed result (see Remark 5.1 in [27]) shows that the 
linear forms on M 3x3 that define the hyperplanes are not canceled by ^4(1), 
because of the —00 entries. 

In [28], we have proved that if the A^s are i.i.d. and the sequence has the 
MLP, then (x(n, X°)) n ^ satisfies the same limit theorem as a sum of i.i.d. 
real variables. Here we prove that it still satisfies the CLT if the A^s are 
mixing quick enough. Quick enough means that the A n 's satisfy the same 
integrability and mixing hypothesis as the real variables in the CLT for the 
sum of stationary variables by Billingsley [6] and Ibragimov [24]. Moreover, 
this proves the CLT under weaker integrability condition than in [28]. 

2.2. Law of large numbers. There have been many papers about the law 
of large numbers for products of random max-plus matrices since its intro- 
duction by Cohen [12]. We can, for instance, cite Baccelli [1], the most recent 
paper by Bousch and Mairesse [8] and Merlet [29] (in French). The latter 
article gives results for a larger class of topical operators, called uniformly 
topical. 

Vincent [31] proved a law of large number for topical operators that will 
do in our case. He noticed that (Vj Xi(n, 0)) n£ N [resp. Xi(n, 0)) n eti] 1S 
subadditive (resp. superadditive), which leads to the following: 

Theorem 2.1 (Vincent [31]). Let (A n ) n& ^ be a stationary ergodic se- 
quence of topical operators and X° an W 1 -valued random variable. If A\.0 
and X° are integrable, then there exist 7 and 7 in M such that 

lim = 7 a.s. and in JL , 

n n 



urn =7 a.s. and in JL . 



1? 



Baccelli and Mairesse give a condition to ensure 7 = 7, hence, the conver- 

neN : 



gence of (-H7 — ' 



Theorem 2.2 (Baccelli and Mairesse [5]). Let (A n ) n£ ^ be a stationary 
ergodic sequence of topical operators and X° an M. d -valued random vari- 
able such that Ai.O and X° are integrable. If there exists an N, such that 
An ■ ■ ■ A\ has a bounded projective image with positive probability, then there 
exists 7 in R such that 

, x(n,X°) , , 

hm = 71 a.s. and in L . 

n xi 
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That being the case, 7 is called the Lyapunov exponent of the sequence. 
Since matrices with rank 1 have a bounded projective image, any ergodic 
sequence (A n ) n ^ with the MLP fulfills the hypotheses of Theorem 2.2. 

2.3. Statements of the results. Let us state the definitions of mixing to 
be used in the sequel. 

Definition 2.2 (Mixing). We denote by T n the c-algebra generated by 
the A^s for k <n and by T n the one generated by the A^s for k>n. We 
define a n and 4> n by the following: 

1. m Q) = S up{ |P(A %y )P(B)l |A Gf.Bee} and n = sup fc cf>(F k , 

2. a(F, G) = sup{|P(ylnS) -F(A)F(B)\\ A eT,BeG} and a n = sup fc a(F h , 

Theorem 2.3. // (A n ) nG fq has the MLP and satisfies one of the follow- 
ing hypotheses: 

A. ^liOGL 2 and Y^=i V<h< +°°, 

B. AiO E h 2+s and Zn=i a n /{2+5) < +00 for some 5>0, 

C. AiO E L°° and £n=i a n < +00, 

-=(x(n,X )-n7l)^JVl, 

where N is a random variable with zero-mean Gaussian law ( or Dirac mea- 
sure in 0) whose variance does not depend on X° , and — > denotes the con- 
vergence in law. 

Moreover, if X° is integrable, then the variance a of M is given by 

1 /2cr 2 \ 1/2 
lim —=E\Jxi(n,X°) — nj = 

re— >+oc Jn v 



7T 

and a = if and only if the sequence (x(n,X°) — wyl) n £^ i- s tight. 

Remark 2.1 (I.i.d. case). When the A n are i.i.d., I gave more precise 
results about a in [28]. In that case, if tf) is a topical function from M rf to 
K, such that sup x \^(A±x) — ip(x)\ has a second moment or if ^4x0 has a 
(6 + e)th moment and X° has a (3 + e)th moment, then: 

• a 2 = lim^E(ip(x(n,X )) -n-/) 2 , 

• a = iff there is a 9 E 7op d with rank 1 such that, for any A in the support 
Sa of Ai and any with rank 1 in the semi-group Ta generated by Sa, 
we have 

6 AO' = 99' + ^1. 



CLT FOR SRS OF TOPICAL OPERATORS 



7 



I also proved that if there is such a 9, then every 9 E T4 with rank 1 has 
this property. 

Moreover, when the A n are defined by matrices in the max-plus algebra, 
a is positive provided that the support of A\ is not included in a union of 
finitely many hyper planes of M dxrf . 

In this paper's framework it is not possible to express a 2 as a limit like 
in the i.i.d. case, because the stationary random variables in Ishitani's proof 
of Theorem 3.1 are not necessarily L 2 (see Section 3.8). 

3. Proofs. 

3.1. Results of Ishitani. We use the results of Ishitani [25] for mixing 
subadditive processes, which we state now: 

Let (fl, T, T, P) be an ergodic measurable dynamical system, and (J r a)a,b&N 
a family of sub cr-algebras of J 7 , such that ^F a t\ = T~ x T h a , and for any 
a < c < d < b, !F™ C T h a . The family (xgt) a <t of random variables is adapted 
if, for any s,t, x s t is ^-measurable. It is subadditive (resp. submultiplica- 
tive) for any s<t<u, x su <x st + x tu (resp. x su < x st ■ x tu ). 

Theorem 3.1 (Ishitani [25] and Hall and Heyde [19]). Assume {x st ) s <t 
is adapted and subadditive. We set T n = J-q and J- n = J-^ 00 " , and define ot n 
and 4> n like in Definition 2.2. We set (p,9) as follows: 




< +00 for some 5 > 0. 



If the following hypotheses are satisfied: 




then 



—^(xon - 717) -^Af, 

where N is a zero-mean Gaussian law (or a Dirac measure in 0). 
Moreover, the variance a of M is given by 
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In the sequel we take = Top^, T the shift and P such that the law of 
(^4n)neN is the image of P by the projection on the positive coordinates. 
From now on, A n is the projection on the nth coordinate, and we denote 
A Q by A, so that A n = A o T n . 

For any s < t, we set x st = ViC^t-i ■ ■■A s 0) i , and T\ = a(A s , . . .,At-i), 
so that (x s t) s ,t£N is adapted to (^2) ,66N- Vincent has noticed in [31] that 
(x s t) s ,teN is subadditive. From now on we check that it satisfies hypothe- 
ses 1-3 with (p, 9) = (2,2) under hypothesis A, (p,6) = (2 + 5, ^rj) under 
hypothesis B and (p,9) = (+oo, 1) under hypothesis C. 

Since x i-> ViG^t-i ' ■ • ^l^Oi is topical, Ait 40 )* 1 < AO < Vi(^0)jl implies 

(2) A( AO )^ x o'- s i^V( A °)»- 

Therefore, we can take VI/ = |^40|oo and hypothesis 2 of Theorem 3.1 is 
checked. In the sequel we check the other two hypotheses. 

3.2. Bound on E(xot) — t-y. It is well known and easy to check that, for 
any A £ Top d and x G M. d , the quantity Vi(^)i — Vi x i onr y depends on A 
and x. We denote it by £(A,x). With this notation, we have 

n-1 

(3) v^(^ x0 )-V( x0 )i=E^'^ x °))- 

i i fc=0 

It follows from the main theorem of [26] — which can be extended without 
difficulty from max-plus to topical operators — that there is a choice Y of 
X°, such that x(n,Y) = Y o T n . In this case, we see that £(Ak,x(k,Y)) = 
£(A,Y)oT k , therefore, V \xi(n,Y) — Vi Y is the partial sum of the stationary 
sequence (£(A, Y) o T fc ) fcgN . 

Let us assume for a while that Y is integrable. Then, so is £(A, Y), because 
AO + Aj Yl < AY <A0 + Vi Yl implies 



V(^o), 



+ |Y 



(4) l^,r)l< 

Therefore, it follows from equation (3) with X° = Y that 

E(\J Xi (n,Y)\ - E (y^) = nE(£(A,F)). 

Since topical functions are nonexpanding, we have [Vi^iC^jY) — x 0n| < 
[[ylloc^herefore, |E(x „) - nE(£(A,Y))\ < 2E(||Y|| 0O ). In that case 7 = 
E(£(.A, Y)) and hypothesis 1 of Theorem 3.1 follows from the integrability 
of Y. 

The end of the subsection is devoted to the proof of the bounds that will 
give this integrability. 
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First, we recall from Mairesse's proof that there is almost surely an n £ N 
such that rk(A_i ■ ■ ■ A_ n ) = 1 and that, for such an n, Y = A_\ ■ ■ ■ A_ n 0. In 
the sequel we denote by N the smallest such n. 

Since (Vj Xi(n, 0)) ne pj [resp. Sj(n, 0)) ne N] is subadditive (resp. super- 
additive), we have, for any n € N and i £ [l,d\, 

n n 
fc=l i fc=l i 

therefore, |^4_i • • -^_ n 0|p < Efc=i |^-feO|p and 

n 

neN k=l 
= E 1 {A r >fc}l" 4 -fcO|p 

= E 1 {r-fc(A_ 1 -A_ fc+1 )^l}l^-fcO|p. 

fceN* 

Finally, we get 

(5) nii<nie<Eii 1 {^(A fc _ 1 ...A l) ^i}i^obL- 

*; 

The finiteness of the right part of this inequality with 1 instead of 6 
would be enough to check hypothesis 1, but the finiteness of this quantity 
also ensures hypothesis 3, as will be shown in the next section. Finally, 
Sections 3.4 to 3.6 will be devoted to the proof of the finiteness under each 
hypothesis of Theorem 2.3. 

3.3. Bound on ||sot — x\t — E(xoj — xi^F^We- We denote by A^ the 
quantity |sot — x\t — E(xoi — #it|A)) • • • >-<4n)|- If t < n, then A^ = 0. From 
now on, we assume t>n. 

First, it follows from equation (2) and the .T^-measurability of AO that 

(6) /\(A0)i<E(x 0t -x lt m i )< V(i0)i 

i i 

and A* < \{A0)\t>. 

Second, if rk{A n ^\ ■ ■ ■ A\) = 1, then 



Xot ~ X U - (X 0n - Sin) = C( A t-l ' ' ' Ai, A n -1 ' ' " A)0) 

(7) 



-£{A t - 1 ---A n ,A n _ 1 ---A 1 0) = 0, 
where £ is the same function as in equation (3). Therefore, we have 

( 8 ) l{rfc(A n _i- J 4i)=l} (xot - Sit) = l{rfc(A„_i-Ai)=l} ( x 0n ~ Sin) 
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and 

l{rfc(A, l _i-Ai)=l} E ( ::c 0t - Zltl-To 1 ) = ^(^{rk(A n - 1 -A 1 )=l}i x 0t ~ Xlt)\Fo) 

(9) = H 1 {rk(A„^ 1 -A 1 )=l}(xOn ~ Xi n )\J^) 

= l{rfc(A„_i-yli)=l}(£0n - Xin). 

Equations (8) and (9) together imply that l{ r k(A n _i— Ai)=i} = 0i an d fi- 
nally, we have 

( 10 ) A n = l{rfc(Ai-i-Ai)^l}^n - 1 {rk(A n _ 1 -A 1 )^l}\(A0)\ V . 

It follows from equations (5) and (10) that (x s t) satisfies hypotheses 1 
and 3 of Theorem 3.1, provided that 

oo 

(H) 5Ill 1 {^(An-i-A 1 )^i}l(^0)|p|| e < OO. 

n=l 

The next three subsections will prove that relation (11) is satisfied, under 
each of the hypotheses of Theorem 2.3. 

3.4. Finiteness under hypothesis A. From the definition of 0, we see 
that, for X G L^JF) and Y G L°°(<7), \E(XY) -E(X)E(F)| < Q)\\X\\ X H^Hoo- 
We apply this inequality with X = |A0|f, and F = l{ r *;(A„-i-A„ /2+1 )^i} ! 
where n/2 is the integer part of the half of n, and we take the square root. 
We get 

||l{rfc(A l -i-"Ai)^l}l^°l'Hl2 

< ll 1 {rfc(A n _ 1 ...^ n/2+1 )^l}l^0|p|| 2 



< ^/F(rk(A n/2 ---A 1 )^l)\\\A0\ P \\ 2 + y/^~ 2 \\\A0\ v \\ 2 . 



The ^(j) n / 2 s are summable by hypothesis A. Let us see that the 

^jF(rk(A n / 2 • • • A\) ^ l)'s are too. For any integers n and no, we have the 
following inequality: 

F(rk(A n+2no -- 

< E (l{rfc(A„-Ai)^l}l{rfc(A, l+2 „ -A n+no+1 )^l}) 

< (<A„ +E(l{rfc(yl n+2no -yl n+no+1 )^l}))E(l{ r .fc(A n -A 1 )^l}) 

< + P(rfc(A no • • • Ai) ^ l))P(rfc(A, • • • Ai) / 1). 

Taking no big enough, we have (4> no + F(rk(A no • • ■ A\) / 1)) < 1, hence, 
F(rk(A n ■ ■ ■ A\) ^ 1) decreases exponentially fast and 



SneN ^(rk(A n/2 ■ ■ ■ A\) ^ 1) < oo. This concludes the proof of hypothe- 
ses 1 and 3 under hypothesis A. 
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3.5. Finiteness under hypothesis B . Let us take X = \A0\^ + ^^ 1+S \ Y = 
1 {rk(A n -A n/2+1 )^i} and q = j^m the mixing inequality (see, e.g., [19]) which 
states for any X G L 1 ^) and Y G L°°(£) 

\E{XY) - E(X)E(Y)| < ^-^{fmx || g ||y||oo 
and let us elevate it to the power i±| . We get 



2+5 

||l{rfc(A„_i-Ai)^l} II (2+<5)/(l+<5) 

< l|l{rfc(A„-i-A n/2+1 )#l} II (2+5)7(1+5) 

< Hrk(A n/2 ■ ■ ■ A{) ± if+V/QW || | A0\ T || (2+5)/ ( 1+ 5) 

+ (6a„ /2 ) 5/(2+5) |||A0|p|| 2+5 . 

The (6a n / 2 ) <5// '' 2+ ' 5 ' > 's are summable by hypothesis B. To see that the 
F(rk(A n/2 ■■■A 1 )^ l)( 1 + <5 )/( 2 +' 5 ) too, we apply the following lemma from [23] 
with A = 2±£ and M st = l{ r k{A t -A s )^i}- 

Lemma 3.1 (Hennion [23]). Let (M st ) s< t be sub-multiplicative and adapted 
th values in [0, 1] sui 
exists c£l, such that 

ln^ A 

n 

This concludes the proof of hypotheses 1 and 3 under hypothesis B. 



with values in [0, 1] such that lim n E(Mo n ) = 0. If *}Z n ull X < oo, then there 



E(M 0n ) < c 



3.6. Finiteness under hypothesis C. We notice that 

X]ll 1 {rMAk-r--Ai)^l}l^°Mli < X^II^CAfe-i-Ai^llilll^Olplloo 
k k 

= \\\A0\ v \\ oo Y,nR>k) 

k 

= \\\A0\ v \\ oo E(R), 

where R = min{re|rfc(^4 n _i • • • A\) = 1}. Moreover, if ¥(rk(A no ■ ■ ■ A\) 7^ 1) < 
1, then i? — no is bounded from above by the hitting time of {rk(A no ■ ■ ■ Ai) = 
1}. The integrability of R will follow from the next theorem due to Chazottes. 



Theorem 3.2 (Chazottes [10]). Let (Q,F,F,T) be a measurable dy- 
namical system, B G J 7 a set with positive probability, and 1b its indicator 
function. If the mixing coefficients a n of the sequence (1b oT n ) ng pj satisfy 
J2n a n < oo> then the hitting time of B is integrable. 
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To apply the theorem, we notice that, when B = {rk(A no ■ ■ ■ A\) = 1}, 
every a n defined by {1b oT n )„, g n is less than the a n -n defined by {A n ) n ^. 
This ensures the hypothesis of Theorem 3.2 and concludes the proof of hy- 
potheses 1 and 3 under hypothesis C. 

3.7. Conclusion of the proof. In the last six subsections we have proved 
that, under hypothesis A, B or C of Theorem 2.3, (Vj Xi(n, 0)) nS N satisfies 
the hypotheses of Theorem 3.1. Therefore, we have 



■M. 



Since topical functions are nonexpanding 



(12) 



-L\/ Xl (n,X°) - Xi (n,0) 
'n v \ n v 



< 



\X l 



n 



0. 



therefore, -^(Vj x i{ n i X°) — wy) — > N. 

o~n = (x(n,X°) — y i Xi(n,X )1) is a function of x(n,X°), which is con- 
verging in law (and even in total variation) by the main theorem of [26], 

therefore, -^S n and A^(x{n, X°) — wyl) — > jVl, which concludes the 
proof of the convergence in law. 
Inequality (12) also implies that 



E 



n 



V^(n,X°)-\/^(n,0) 



< 



E||X L 



0, 



so that 



1 



lim — E 



\Jxi(n,X°)-nj 

i 

follows from lim„^ +00 -^E|x 0n - ny\ = (^) 1//2 



TT 



3.8. Tightness. Without loss of generality, we assume 7 = 0. (Otherwise, 
just replace A n by A n — 7.) One part of the equivalence is obvious. To 
prove the other part, we have to go into the proof of Theorem 3.1. Ishitani 
constructs a random variable Z (named yoi m [25]) and approximates XQ n 
by the Birkhof sum 

Sn = TJl=o Z ° Tk (Von in [25]). Then he shows that 
(5 n ) ng N fulfills the hypotheses of Billingsley-Ibragimov's CLT. 

In Billingsley-Ibragimov's CLT, the asymptotic variance is zero if and 
only if Z is a coboundary, that is, if there is a random variable / such that 
Z = foT-f. (See, e.g., [19].) 
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Let us assume we are in this situation and identify Z. It is built as a kind 
of Cesaro type limit of the sequence (xon — £in)neN- But in our situation 
equation (7) says that this sequence is ultimately constant and that X n is 
equal to the limit as soon as rk(A n ■ • ■ A\) = 1. 

Let us denote by R the smallest such n and by ip the topical function that 
maps x G l rf to V/i x i- The random variable R is almost surely finite because 
of ergodicity and MLP. With notation, we have 

Z = x r - x m = tp(A R ■ ■ ■ A 0) - ip(A R ■ ■ ■ 40) a.s. 

and, for any integer n such that rk{A n • • -A\) = 1, 

(13) foT-f = ^(A n ---A 0)-iP(A n ---A 1 0). 

In the sequel we deduce the tightness from equation (13). As a first and 
main step, let us show that {^{Ar ■ ■ ■ A_ n 0)) n ^ is tight. For any k G N, 
since rk(A R ■ ■ ■ 4) = 1, rk(A R ■ ■ ■ A_ k ) = 1 and equation (13) holds for n = 
R o T k + k. Compounded by T~ k , it becomes 

/ o T- k+l -foT~ k = ^{A R ■ ■ ■ A^ k 0) - ^{A R ■ ■ ■ A_ k+1 0). 

Summing over k, we get 

/ o T - f o T~ n = ^(A R ■ ■ • A_„0) - ^(A R • • • AoO), 

from which the tightness of (^P(Ar ■ ■ ■ A_„0)) ng N is obvious. 

The tightness of (A_i ■ ■ • ^4_ n 0) ng N is obvious too, because the sequence 
converges in law. 

From those two tightnesses, we successively deduce the tightness of the 
following sequences: 

- (ip(A_i ■ ■ -^4_ n 0)) ne N, because equation (4) implies that 

\i>{A R ■ ■ ■ A_ n 0) - V(A_x • • • A_„0)| = \£(A R ■ ■ ■ 4,,A_i---A_ n 0)| 

< \i/>(A R ■ ■ ■ 4)0)| + \A- 1 ---A. n Q\ v . 

- ((ifj(A-x ■ ■ ■ A_ n 0, A-i ■ ■ ■ A_ n 0)))„ e N, again because (A_i ■ ■ ■ ^_ n 0) n6N is 
tight. 

- (A-i ■ ■ ■ ^4_ n 0) n gN, because x \— > (ip(x),x) is a bi-Lipschitz homeomorphism 
from M. d to R x PM^ ax (see [28]). 

- (x(n, 0)) ne N 5 because, for any n £ N, the random variables (^4_i • • • A_ n 0) 
and x(n,0) have the same law. 

- Eventually (x(n,X°)) n£ ^, because the A n are nonexpanding and, there- 
fore, we have \\x(n,X°) — x(n,0)||oo < ||^°||oo- 

Acknowledgments. This work was done while I was a JSPS post-doctoral 
fellow at Keio university. I gratefully thanks H. Ishitani for introducing me 
to the article [25] and for interesting discussions, and M. Keane for his useful 
rereading. 



14 



G. MERLET 



REFERENCES 

[1] Baccelli, F. (1992). Ergodic theory of stochastic Petri networks. Ann. Probab. 20 
375-396. MR1143426 

[2] Baccelli, P., Cohen, G., Olsder, G. J. and Quadrat, J. P. (1992). Synchroni- 
sation and Linearity. Wiley, New York. MR1204266 

[3] Baccelli, F. and Hong, D. (2000). Analytic expansions of max-plus Lyapunov 
exponents. Ann. Appl. Probab. 10 779-827. MR1789980 

[4] Baccelli, F. and Hong, D. (2000). Analyticity of iterates of random non-expansive 
maps. Adv. in Appl. Probab. 32 193-220. MR1765163 

[5] Baccelli, F. and Mairesse, J. (1998). Ergodic theorems for stochastic operators 
and discrete event networks. In Idempotency {Bristol, 1994) 171-208. Cam- 
bridge Univ. Press. MR1608402 

[6] Billingsley, P. (1968). Convergence of Probability Measures. Wiley, New York. 
MR0233396 

[7] Blondel, V. D., Gaubert, S. and Tsitsiklis, J. N. (2000). Approximating the 
spectral radius of sets of matrices in the max-algebra is NP-hard. IEEE Trans. 
Automat. Control 45 1762-1765. MR1792275 
[8] BOUSCH, T. and Mairesse, J. (2006). Finite-range topical functions and uniformly 

topical functions. Dyn. Syst. 21 73-114. MR2200765 
[9] Braker, H. (1993). Algorithms and applications in timed discrete event systems. 
Ph.D. thesis, Delft Univ. Technology. 

[10] Chazottes, J.-R. (2003). Hitting and returning to non-rare events in mixing dy- 
namical systems. Nonlinearity 16 1017-1034. MR1975794 

[11] Cohen, G., Dubois, D., Quadrat, J. P. and Viot, M. (1985). A linear system 
theoretic view of discrete event processes and its use for performance evaluation 
in manufacturing. IEEE Trans. Automat. Control AC-30 210-220. MR0778424 

[12] Cohen, J. E. (1988). Subadditivity, generalized products of random matrices and 
operations research. SI AM Rev. 30 69-86. MR0931278 

[13] Crandall, M. G. and Tartar, L. (1980). Some relations between nonexpansive and 
order preserving mappings. Proc. Amer. Math. Soc. 78 385-390. MR0553381 

[14] Gaubert, S. and Hong, D. (2000). Series expansions of Lyapunov exponents and 
forgetful monoids. Technical report, INRIA. 

[15] Gaubert, S. and Mairesse, J. (1999). Modeling and analysis of timed Petri nets 
using heaps of pieces. IEEE Trans. Automat. Control 44 683-697. MR1684424 

[16] Gaubert, S. and Mairesse, J. (1998). Task resource models and (max, +) automata. 

In Idempotency (Bristol, 1994) 133-144. Cambridge Univ. Press. MR1608394 

[17] Gaujal, B. and Jean-Marie, A. (1998). Computational issues in recursive stochas- 
tic systems. In Idempotency (Bristol, 1994) 209-230. Cambridge Univ. Press. 
MR1608327 

[18] Gunawardena, J. and Keane, M. (1995). The existence of cycle times for some non- 
expansive maps. Technical Report HPL-BRIMS-95-003, Hewlett-Packard Labs. 

[19] Hall, P. and Heyde, C. C. (1980). Martingale Limit Theory and Its Application. 
Academic Press Inc, New York. MR0624435 

[20] Heidergott, B. (2000). A characterisation of (max, +)- linear queueing systems. 
Queuemg Syst. Theory Appl. 35 237-262. MR1782609 

[21] Heidergott, B. and de Vries, R. (2001). Towards a (Max, +) control theory 
for public transportation networks. Discrete Event Dyn. Syst. 11 371-398. 
MR1852744 

[22] Heidergott, B., Oldser, G. J. and van der Woude, J. (2006). Max Plus at 
Work. Princeton Univ. Press. MR2188299 



CLT FOR SRS OF TOPICAL OPERATORS 



15 



[23] Hennion, H. (1997). Limit theorems for products of positive random matrices. Ann. 

Probab. 25 1545-1587. MR1487428 
[24] Ibragimov, I. A. and Linnik, Yu. V. (1971). Independent and Stationary Sequences 

of Random Variables. Wolters-Noordhoff Publishing, Groningen. MR0322926 
[25] Ishitani, H. (1976/77). A central limit theorem for the subadditive process and 

its application to products of random matrices. Publ. Res. Inst. Math. Sci. 12 

565-575. MR0438475 

[26] Mairesse, J. (1997). Products of irreducible random matrices in the (max, +) alge- 
bra. Adv. in Appl. Probab. 29 444-477. MR1450939 

[27] Merlet, G. (2004). Memory loss property for products of random ma- 
trices in the (max, +) algebra. Technical report, IRMAR. Available at 
http : //hal . ccsd . cnrs . f r/ccsd-00001607 . 

[28] Merlet, G. (2005). Limit theorems for iterated random topical operators. Technical 
report, IRMAR. Available at http://hal.ccsd.cnrs.fr/ccsd-00004594. 

[29] Merlet, G. (2005). Produits de matrices aleatoires: Exposants de Lyapunov pour 
des matrices aleatoires suivant une mesure de Gibbs, theoremes limites pour 
des produits au sens max-plus. Ph.D. thesis, Univ. Rennes. Available at 
http : //tel . ccsd . cnrs . f r/tel-00010813 . 

[30] Resing, J. A. C., de Vries, R. E., Hooghiemstra, G., Keane, M. S. and Ols- 
der, G. J. (1990). Asymptotic behavior of random discrete event systems. 
Stochastic Process. Appl. 36 195-216. MR1084975 

[31] Vincent, J.-M. (1997). Some ergodic results on stochastic iterative discrete events 
systems. Discrete Event Dynamic Systems 7 209-232. 

CEREMADE 

Universite Paris-Dauphine 

Place du Marechal de Lattre de Tassigny 

75775 Paris cedex 16 

France 

E-MAIL: glenn.merlet(Q!gmail.com 



